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QJ I Abstract 

^ ■ In [1], a two level model of the occurrence of financial pyramid (bubbles) has 

, been considered. We also considered the mathematical analogy of this model to 

Bose condensation. In the present paper, we explain why Ponzi schemes and bubbles 
result in a crisis in real economics. In [2j , the law of increase of entropy in financial 
systems, and consequently increase of Kolmogorov complexity, is formulated. If this 
law is broken, the financial system makes a phase transition to a different state. In 
^ ' [3] the author considered a two level model of the zeroth-order phase transition 

^4H . which was interpreted in [Jj as an analog of social catastrophe. In the present paper 

Q^l we also examine this model. 



By the end of 2008, the United States has accumulated considerable debts. The debts 
! can be counted by the same scheme as GDP, that is, as the difference of non-repayed 
^ \ and repayed debts in a year provided that this difference is not zero, i.e., the debt has 



■ increased owing to the velocity of debt. 

Just as negative money, i.e., debts, do not occur in the definition of velocity of money 
and in the calculation of the GDP, so the definition of the velocity of debt does not involve 
OO . positive money (earned money, etc.). 

^ I If Ni is the debt resulting from one turnover (one act of borrowing from somebody 
and paying the dept to somebody else), then the total debt is equal to Mj = Nia, where 

^ ■ cr, the number of turnovers per annum, will be called the velocity of debt (by analogy 
. with "velocity of money"). The highest velocity of debt is observed for bubbles, or Ponzi 
schemes. A mortgage debt has only one turnover. For a car loan, the velocity of debt is 
also equal to one if the customer made a large borrowing and has not repaid the loan yet. 

Thus, 'Y^NiCTi < M, where M is the total debt (like a "negative GDP"). Hence the 
debts obey the same theorems of number theory as GDP (see [1] and also [5]). 

If M is the total debt and a is the average velocity of debt (analog of the temperature 9 
in physics), then the "negative money supply" is given by iV = M/a, where N = Y. .=i 
/c ^ 1, and Ni corresponds to the ithe set of subjects, which has the cardinality a,. If 
we adopt, similarly to the Baumol-Tobin model, that the velocity of debt of a single 
subject is proportional to the square root of its debt (cf. [TOl Appendix to Chapter 13]) 
and take into account that the debt is equal to iNi, then is proportional to i and it 
is more accurate to assume that Ni = irii, where rij are distributed uniformly. Hence 
Yl'i=i Oii'^'^f^i = M and ^ irii = N. For we assume that they decrease according to the 
Pareto distribution. 

The critical number Nq after which the "Bose condensation phenomenon" occurs is 
computed in [H [TTl [M]. In these papers, the Pareto distribution was proved considered 
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as negative dimension, namely, it was proved that the larger the Pareto index for the 
distribution of and the lesser the Baumol-Tobin dependence, then the larger the critical 
number Nq and the more difficult it can be overcome. 

However, there is a substantial difference between inflation occurring if N > Nq for 
the case of income and the case in which N > Nq for debts. 

In the first case (infiation) small denominations fall out of circulation (part of the 
money is abolished), and the nominal value of money is reduced. In the second case, 
debts cannot be annihilated by denomination, and excessive debts for N > Nq accumulate 
at the "lower level", i.e., at the subjects who have one turnover per annum. The zero 
velocity of debt will be taken into account in incomes rather than debts. On the other 
hand, it follows from the theorem proved in [TOl Appendix to Chapter 13] that small 
velocities of both the money and the debts practically do not affect the value of Nq. 
Hence the definition of nonzero lower level of velocity of money is not essential in the 
rather schematic GDP model. 

However, since the lower velocity level for debts plays an exceptional role and the 
model involving the number of turnovers is too coarse for determining the least velocity 
of debt, one has to single out the "longest-term debts" whose repayment spans over many 
years, i.e., whose repayment velocity is small. 

How the long-term project contribution to GDP is calculated? It is necessary to 
calculates the excess of the cost of the project already performed part in a year over the 
cost of the entire project in L years minus the total expenditures in L years divided by L. 

Similarly, the rate of the debt "turnover" in the case of a long-term credit obtained 
for L years is equal to 1/L. Therefore, if, without loss of generality, we assume that the 
numbers of turnovers larger than 2 are integer numbers, then, for a single "turnover" (a 
credit taken once), we must take the credit period (the number of years L) into account. 
Correspondingly, the rate of the debt turnover equal to 1 splits into rates of the debt 
turnover that are less than 1. Hence the debt Bose condensation accumulates precisely at 
these slow rates of the debt turnovers, and the bankruptcy expects precisely them. This 
follows from the mathematical theorem of the number theory. But the theorem does not 
speak how this is reahzed in the society. Possibly, to escape the bankruptcy, the debtors 
with high level of the debt velocity must again take long-term deposits. Perhaps, this 
chain can be explained by economists, but the mathematical law inexorably predicts this 
scenario of the debt crisis. 

It is on these debts that the excessive debts exceeding Nq mainly fall; as a result, the 
long-term borrowers default on their debts, which fall out of the chain and descend on 
the subsequent levels. 

Hence mortgage is the first to collapse, and the car loans follow. 

Thus, the mathematical theory, which pertains to number theory, is the same, but the 
result in economics is different. 

It follows from the preceding that to fight the crisis and to avert it in future one 
should favor a debt structure in which the share of high-velocity debts is smaller that 
in the existing structure. In [H |TT], it was proved that if the Pareto index is not too 
large, then it is possible to increase the critical number by introducing a larger number 
of currencies, for example, in each of the United States, by introducing its own currency 
in addition to the dollar. 
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1 Economic Law of Increase of Kolmogorov Com- 
plexity 



In his celebrated paper on complexity, A.N. Kolmogorov defined Kolmogorov complexity 
as a minimal code of a sequence of numbers. He also mentioned some other complexity, 
which was related to decoding. For the stock exchange, the definition of a code length 
that involves decoding is more appropriate to the situation under consideration, which 
leads to an arbitration-free situation in the end. 

Suppose that a trader contrived a combination that brought a big gain to this trader. 
The market intuitively deciphers this combination (algorithm) and, thereby, annihilates 
it. Then, the trader invents some more complicated algorithm, and so on. In the end, the 
complexity of the algorithm becomes so high that it only slightly differs from a pseudo- 
random complexity. In the limit, this leads to an arbitration-free situation. 

For this reason, it is more natural to understand complexity as the length of a code 
with polynomial decoding complexity. This code should be used to encode a sequence. 
The longer the code, the closer the situation under consideration to general position, 
which corresponds to the law of large numbers. If the minimum specified above tends to 
zero as ^ oo, then the set of deals {Ni} is close to an arbitration-free set. 

In the language of probability theory, the initial probabilities corresponding to Aj, 
which equal Gi/G in the sense of von Mises, tend to zero as — > oo (and, hence, as 
G^oo). 

Therefore, the initial mathematical expectation 

T.9^ ~ ^ G 

is such that each term of the probability vanishes in the limit. However, the mathe- 
matical expectation tends to a constant larger than zero and the partial sums Ei=i Vi 
do not tend to zero. 

Similarly, the mathematical expectation ^ does not tend to zero, whereas ^ — 0, 
but 4? y^f" • Aj does not tend to zero. 

This is a very important economy factor, especially in the psychological aspect (see 

glE]). 

Since the tendency of the market to an arbitration-free situation is equivalent to the 
eagerness of the market participants to gain profits, it follows that the entropy, if it 
depends on some parameters, must increase with respect to these parameters most rapidly, 
in the direction of steepest ascent. 

If li{E, N) depends on k additional xi, ■ ■ ■ , Xk, i.e., H = H{xi, ■ ■ ■ , Xk), (the depen- 
dence on E and A^ is not indicated), then 



'dt 



For simplicity, suppose that C{H) = 1. Otherwise, setting F{H) = f C{H)dH, we 

dt 



obtain ^ = (V-F)^. This implies 



(x - 

H{x, t) = min ( ^ ^ + H{i) 



i V 2t 



Remark 1. Performing geometric quantization [7], we obtain 



(x,t) = ln-i=|e-^e^(«)rfe, 



This asymptotically coincides with the tunnel canonical operator on the Lagrangian man- 
ifold Pi = ^ shifted along the trajectories x= p and P= 0. The quantization of Eq. ([T]) 
leads to solution of the heat equation [S] 



Quantizing this equation (as a Bose system) and setting u = \1'^ and = \& ■ \E'~, we 
obtain a linear equation whose asymptotics is again the canonical operator. 
The prices Aj(xi, ■ ■ ■ ,Xk) vary according to the equation 



Thus, if we know the change of some (specific) price, we can find the function C{H) 
and, thereby, solve the problem completely. 

A similar situation may arise when, e.g., barriers between countries (such as customs 
tariffs) decrease and prices tend to even according to the law of steepest ascent of entropy, 
because profiteers begin to use the arising arbitration (difference in prices) at a furious 
pace, which leads to an arbitration-free situation most rapidly. Such was the situation in 
the Soviet Union in the late 1980s when the iron curtain collapsed 

On September 30, 2008, the National Debt Clock in Manhattan ran out of digits as 
the United States public debt exceeded $10 trillion, a significant symbol of the current 
financial and economical crisis. There are many diverse factors behind the crisis, and 
they are naturally a subject of broad interest. However, apart from economics laws, 
which are widely discussed nowadays, there is a purely mathematical law contributing to 
the disaster, an inexorable law of numbers, which economists, let alone the general public, 
fail to recognize. 

Suppose we want to deposit two kopecks in two different banks (see [1]; then we can 
say that there are three possibilities: (1) put both kopecks into one of the banks; (2) put 
both kopecks into the other bank; (3) one kopeck in one bank, the other, in the other 
one. Here it is of no consequence which of the two coins we deposit in the first bank and 
what is its year of issue. Now imagine a situation in which we are depositing one kopeck 
and one pence instead of two kopecks. In that case, we have four options rather than 
three, because it is significant which coin we placed in what bank, and so the variant in 
which the coins are placed in different banks yields two different options: (a) one kopeck 
to bank 1 and one pence to bank 2; (b) one kopeck to bank 2 and one pence to bank 1. 

These two cases give rise to two different laws in number theory (see P]). In the 
second case, the corresponding law implies that the total number of cents (and hence of 
dollars and of millions of dollars), i.e., the total amount of money in circulation, cannot 
be arbitrarily large. There is a threshold, and if the amount of money in a given currency 
exceeds that threshold, then an economic disaster occurs. The same thing pertains to 
negative amounts of money, i.e., debts. 

Mathematically, this kind of disaster can be described as an analog of "Bose conden- 
sation" : should there be too many particles in a system, all the excessive particles would 
collapse into the ground state. The condensation effect is indeed observed in finance: for 





(2) 




(3) 



example, if the inflation rate is too high (too much money has been issued), then the 
lower denominations (like cents, pence, or kopecks) die away, i.e., are withdrawn from 
circulation. 

It is important to note that the condensation phenomenon is solely due to the fact 
that particles are indistinguishable — there is no condensation at all if the particles are 
distinguishable (and so obey the Gibbs statistics). One should simultaneously use as 
many currencies as possible. For example, if, in addition to dollar, its own currency were 
introduced in each of the United States, then the sum of the thresholds for all states would 
be much greater than the threshold in the case of a single currency; it would be much 
harder to exceed this new threshold, and hence the disaster would at least be postponed. 
By the same pattern, all national currencies should have been retained along with the 
euro. Then the violation of a country's threshold could cause a crisis in that country but 
would not affect the other economies, just as the 1998 default in Russia did not bring 
down the world economy. (Unfortunately, the dollar largely circulates outside the United 
States, and hence the world on the whole is much more affected this time). 

Thus, the introduction of K currencies has raised N^r by the factor of approximately 
K^^"^. Or, the introduction of concurrently existing K currencies increases the crisis 
threshold by the factor K^^"^. The outcome will be more or less the same regardless 
of whether these new currencies are nationwide or function only within their respective 
states. 

The above described transition from the MMM pyramid ("bubbles") leads to entropy 
increase, therefore, to the Kolmogorov complexity. Hence "bubbling" leads to decrease of 
complexity. Likewise, transition to only one currency leads to decrease of complexity, and 
this is contributing to a crisis. Development of advanced technologies results in increase 
of complexity, therefore it hinders a crisis. 

Thus, a general law of economics is the law of increase of complexit-^ Otherwise, a 
crisis is imminent. 



2 Nonstationary financial averaging 

Consider a random variable A which takes / different values Ai, . . . , A^. In what follows, we 
assume that the values of this random variable satisfy the following condition: if fci, . . . , /cj 
are integers such that 

I I 
^ A;, = 0, ^ \ki = 0. (4) 

i=l 1=1 

then /cj = for all z = 1, 2, /. 

In |13], the notion of financial averaging and financial averaging axioms were intro- 
duced. Suppose given a convex function f{x) having inverse f~^{x), where x G R and a 
set of nonnegative numbers Pi, where z = 1, we call them weights in what follows. 

The first axiom [13] is that the financial averaging of the random variable A corre- 
sponding to the function f{x) and the set of weightspj, i = 1, . . . , / is the expression 

E/,p(A) = rMX]p./(A.)]. (5) 



The analogy to the above-mentioned effect in statistical physics is given e.g., in [12]. 



We modify the fourth axiom as follows: For any set of weights pi, any set of numbers 
Xi , and any number a, the following equality holds: 

E;,,(A) = Ca + E;,,(A), (6) 

where Aj = Aj + a, and C is a number not depending on a. Then, the following assertion 
is valid. 

Proposition 1. Condition ^ implies that the function f{x) has the form 

f{x) = Aexp{-(3x), or f{x)=Ax + D, (7) 
where the numbers (3 and A are nonzero. 



This assertion is proved by the method described in p^. 

Hereafter, we consider the case of a nonlinear function f{x). According to Propositon 
1, the financial averaging of a random variable A corresponding to a set of weights Pi has 
the form 

E,(A) = -iln|^gp,e-^^^j. (8) 

Properly, Axioms 2 and 3 [T3] do not hold; thus, we regard this averaging as non- 
stationary and consider the process stabilizing it. The stabilization procedure does not 
essentially use the fourth axiom. 

Following [13], we associate the random variable under consideration with an n- 
dimensional space L and an operator A in this space. The operator A has eigenvectors 
Gi, where z = l,2,...,Z, which form an orthonormal basis in the space L. The eigenvalue 
corresponding to the eigenvector ei equals Aj. Consider a series of M identical trials whose 
outcomes are occurrences of the values Ai, . . . , A^ of the random variable. According to 
[H], this probability problem corresponds to the space Cm being the tensor product of 
M copies of L. The vectors in the space Cm have the form 

I I 

^ = ^ ■ ■ ■ ^ ip{ii, iM)ei, ® ■ ■ ■ ® Bi^,^, (9) 

«1=1 «A/ = 1 



where ip{ii, . . . , im) is an arbitrary function of the set of discrete variables ig = 1, . . . , /, 
s = 1,...,M. 

To the series of trials under consideration, there corresponds, in addition to the space 
Cm, the operator Hm acting on vectors ([9]) as 

I I / M \ 

Um^ = ■ ■ ■ ( Z J ^^'1' • • • ' ® ■ ■ ■ ® ^iM- (10) 

n=l iA/=l \s=l / 
In further considerations, we use the notation 

M 

Mi{ii, . . . ,iM) = ^Sii^, (11) 

where da/ is the Kronecker symbol. Let us define projectors P{m} 
I I / I \ 

-P{A/}* = Z ■ ■ ■ Z I n ^M,Mdh,...Ml) j ^"(^1, ■ ■ ■ , ^M)e^, ® ■ ■ ■ O e^^. (12) 

«l=l «M=1 \i=l / 



Here, {M} denotes an arbitrary set of nonnegative integers Mi,...,M„ satisfying the 
condition 

I 

^Mi = M. (13) 

i=l 

It is easy to show that operators f|T3|l satisfy the equahty 

P{m}P{m'} = W5m,m'P{m} (14) 



i=l 



thus, these operators are indeed projectors. The operator TIm defined by ( fTOll is expressed 
in terms of the projectors P{m} as follows: 



{M} \i=l J 



'{MY (15) 



In the subspace >C{m}, onto which the space Cm is projected by P{m}i "we distinguish the 
basis formed by vectors of the form 

En,...,i,,, =ei, (g)---®ei^,, (16) 

where occurs Mi times, or, in other words, the set ii, . . . ,iM satisfies the conditions 

Mi{t^,...,tM) = Mi. (17) 

Consider the following system of vectors ${a/}: 

I 



${A/} = ^{A/}e ® • ■ ■ ® e, e = J2^i- (18) 

i=l 

Vectors ffTSj) are equal to the sum of all vectors ffTB]) from the subspace C^m}] in addition, 

*W = ^ f ri^^' ) Symm,^,...,,^E,,...,^, (19) 



where Symm^^ is the operator of symmetrization over the indices ii, . . . ,iM- 
We endow the space Cm with the norm 

I I 

ii^ii = E--- E m^i,---,^M)\. (20) 

3 The evolution process 

In this section, we define the notions of evolution and data reduction. 

Definition 1. The one-step evolution of an ensemble consisting of M elements, i.e., of a 
vector E Cm, is the following vector belonging to the space Cm- 



^1 = exp ( -/3HAf ) ^, (21) 

where /? is a parameter. 
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Definition 2. The data reduction, or factorization, is tlie operation tliat maps eacli vector 
\I/ G Cm to anotlier vector according to tlie rule 

R{^) = Y,\\P{M}n^{M}. (22) 

{M} 

wliere summation is over all sets of nonnegative integers Mi, . . . ,M„, satisfying condi- 
tion ([H. 



To each set of nonnegative numbers where i = 1, . . . ,1 we assign the following 
vector of the space Cm- 

^3 = g ® • ■ ■ ® g, (23) 
where g is the vector from L, defined by 



I 



g = $^^7^ei. (24) 



i=l 

Proposition 2. If gi = pie~^^\ then the expression 

1 



M(3 '-'"'-'^ 

coincides with financial averaging ^ 
Proof. Let us write vector fl23l) in the form 



In (11^,11) (25) 



^9= [Yl ® ■ ■ ■ ® 5Z s^M^^M = 5Z ■ ■ ■ 5Z • • • 3iM^^^ ® ■ ■ ■ ® Gi^i- (26) 



Substitution of ([2SD into ([201) yields 



M 



\\'^9\\=T.■■■T.9n■■■9^M=[J29^) ■ (27) 
Taking into account that g^ == pie~^^^ and applying fl271) . we obtain 

-ilnf^Ke-"*-). (28) 



The right-hand side of (!28|) coincides with the right-hand side of ([8]), which proves the 
statement. □ 

In what follows, we do not use the axiomatics of nonlinear averaging; instead, we 
consider the evolution of an arbitrary (nonnormalized) probability distribution of gi. 

For any set of nonnegative g^, where i = 1, . . . , / we define the evolution of vectors 
\l/g(n) G C^M} by the recursive formulas 

^,(n+l) = /?(exp(-/57^A/)^p(Ti)), n = 0,l,..., ^,(0) = ^,. (29) 
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We use ipgi^ii, . . . ,iM',n) to denote the coefficients in expansion ([9]) of the vector '^g{n) in 
basis ( fT6l) . Thus, ipgih, ■ ■ ■ ,iM',n) and "^gin). are related by 



I I 

^sH = X] ■ ■ ■ XI • • ■,iM;n)ei^ ® ■ ■ ■ ® eij^^. (30) 

The function ipglii, . . . , ^M; '^) can be regarded as a nonnormahzed distribution function 
of the factored ensemble after n steps under the condition that, at the initial moment, the 
distribution function of the factored ensemble has the form ipgi^ii, . . . , iM', 0) = gi-^ . . . gi,^j. 
In this case, gi has the meaning of a nonnormahzed distribution function for one system 
in the factored ensemble. 
Consider the functions 

F(n,^,M) = - ^^^^^^^ ln(||vl/,(n)||), (31) 

and 

Wi{n, g, M) = 5Z ■ ■ ■ 5Z ^9(^5 ^2, • • • , ihi; n). (32) 

Since ipgiii, ■ ■ ■ iIm) is a distribution function for the factored ensemble of M systems, 
formula ( |32l) defines a distribution function of one system. Consider the limit of distribu- 
tion functions (13T]) and ( |32l) for one system of the factored ensemble as the number M of 
its elements tends to infinity. 

Theorem. As M ^ oc the functions F{n, g, M) and Wi{n, g, M) approach the following 
limits: : 

Jim^F(n,^,M) = -iln (^^e-^^^^"^-^') g]'^-^^)^ ^ F{n,g), (33) 

g-n/3A,/(n+l)^V{«+l) 

lim Wi{n, g, M) = — ' = w,{n, g). (34) 

Proof. Vector ( l23l) can be represented in the form 



{M} 

Substituting (l35i) into ( l29l) . we obtain 



|M} \j=l / 



(35) 



This implies 

n (Mj)» )' ^^^^ 



where are expressed in terms of ii, ... ,iM according to ffTTl) . Substituting fl37j) into fl20|) 
and summing ( 1371) over all but one subscripts, we obtain 

II^.HII = (^0""^En\pv^' (38) 

{A/} i=l ^ ''^ 



42 = 1 iM = l {M} \i=l ^ ' / \j = l 



To determine the asymptotic behavior of the sums in (138|) . (139|) as M — oo, we apply the 
Laplace method. This method gives the equalities 

^hn^l:ln(l|vl/,,,H||) = (n + l)ln ( ^ e-"^^'/("+i)^,^/("+^) ) , (40) 



,i=l 

' -n/3A,/(n+l) 



lim T— -— -- ■ ■ ■ "ippAh "^2, • • • , in; n) = — ; — 77^-777- (41) 

i^' 12 = 1 lM = l Z^]=l^ iJj 

The substitution of ^ and (SI]) into (EI]) and ([32]) yields ([33]) and (JSl]). This completes 
the proof of the theorem. □ 

Thus, by virtue of flMl) . the normalized distribution Wi[0,g) = transforms into 

Wi{n,g) in n evolution steps. 

Now, consider the behavior of F{n,g) (l33l) and Wi{n,g) ( IMl) as n — >■ 00, i.e., the limit 
distribution. Let / be a nonempty subset of the set 1, . . . , /. Suppose that the gi satisfy 
the condition 

gi > 0, for i e I, g, = 0, for i ^ I. (42) 
According to fl33l) and flMl) . for any g^ satisfying fjl2l) . we have 



lim F(n,(7) = -^ln( 



(43) 



lim Wi{n,g) = —, — = pi. (44) 



-/3A 



If the set / coincides with the set of 1, . . . , /, expressions f l^ and (jS]) are, respectively, 
the free energy and Gibbs distribution at temperature of a system with nondegenerate 
energy levels Ai, . . . , A;. If / is smaller than 1, . . . , /, then expressions ( 143!) and (jH]) describe 
the distribution of the system only over a part of levels, rather than over all levels. Note 
also that, if the gi satisfy the condition 

f Ae~^^' for i e / 
= I for . ^ / ' (^^) 

where A is an arbitrary positive number, then Wi{n,g) flMj) coincides with pi (jH]) for 
any n. Thus, the full and partial Gibbs distributions are invariant with respect to the 
evolution described above. 



in 



Remark 2. Consider the special case where the random variable A takes the same values 
as the polynomial 

£{N,)=J2a,NI (46) 

q=0 

where A^i ranges over 0, 1, . . . , A^. Consider the satisfiability of condition (jl]) in this case. 

Proposition 3. Ifp < N, then condition ^ holds for no values of the coefficients Aq, and 
if p > N, then condition holds for generic coefficients Aq. 



n 



4 Zeroth-order phase transition 



4.1 Zeroth-order phase transition: general considerations 

4.1.1 Main large parameters in thermodynamics 

Thermodynamics studies steady-state processes in which, independently of its initial state, 
the system comes to a state that remains the same in what follows. But if such a state 
still varies, then this is a thermodynamic variation only if the process occurs extremely 
slowly [15]. In other words, we slightly vary the state of the system and wait until it 
returns to equilibrium and again is independent of the initial state. 

Thus, a large parameter, time, is invisibly present in this case, i.e., we observe the 
system at huge time intervals. On the other hand, thermodynamics is the limit macro- 
scopic theory obtained from the microscopic statistical physics as the number of particles 
N tends to infinity. Hence there are two large parameters, and many things depend on 
their ratio. In turn, classical statistical physics is the limit of quantum statistics as h ^ 
{h is the Planck parameter). Thus, already three large parameters "meet" here: time 
scale, number of particles, and 1/h. 

In thermodynamics, sufficiently large time periods are also often considered, and in 
these periods only a part of the system comes to equilibrium; this is the so-called "local 
equilibrium." For example, when plasma is heated by a magnetic field, the ions begin 
to obey the Maxwell distribution after a large observation period, and the entire system 
comes to the thermodynamic equilibrium only after a significantly large period of time. 

4.1.2 The Gibbs formula 

Let us consider the Gibbs formula, which has the following form in quantum statistics. 

Suppose that a system is characterized by the Hamiltonian operator if in a Hilbert 
space (in particular, in the Fock space) Suppose that the operator H has a nonnegative 
discrete spectrum Aq, Ai, . . . , A„ . . . . Then the free energy depending on the temperature 9 
is determined as 

oo 

E = e\nJ2e'^-/'^6n, 

n=0 

where Sn is the multiplicity of the eigenvalue A„. 

We note that one more parameter appears in the Gibbs postulate, namely, 

M 

E = e\n lim Ve-^"/^5„. 

M^oo ^-^ 
n=0 

This is one of the most important facts, since limM^oo and limTv^oo do not commute! It 
turns out that one must first pass to the limit as M — oo and then as ^ oo. 

If we speak about local equilibrium, then, as a rule, the Gibbs formula deals with 
a subset {A„'} of the set {A„}. This must be done in the study of the above problem 
concerning the ion distribution in plasma on heating. This is also done in glass and in 
several other physical problems. 

4.1.3 "Friction" operator 

At temperature 6* = 0, the Gibbs formula gives the value of the lowest eigenvalue, which, 
in physics, is called the "ground state." 
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This is a law which has not been stated accurately and which was called the "energetic 
efficiency" by N. N. Bogolyubov and others. The lower energy is energetically more 
efficient. For example, if a system is perturbed by an operator \^ ($ — > $) and this 
operator is sufficiently small, then the transition matrix from the state A„ to the state 
is determined by the matrix element (\E'„V^\i/^), where and are the eigenfunctions 
corresponding to A„ and A^, respectively (see [H]). 

The square of this matrix element is the probability of transition from the state with 
energy A„ to the state with energy Am. But if A„ < Am, then this transition is "not 
energetically efficient" and hence, as if from the viewpoint of statistical physics and ther- 
modynamics, is unrealizable, i.e., from the mathematical viewpoint, it must be set to 
be zero. Then the transition matrix is not a self-adjoint matrix with zero entries above 
the diagonal. This means that friction is taken into account. A pendulum oscillates and 
stops in the end if friction is taken into account, i.e., it comes to the "ground state." The 
"friction" operator was considered in more detail in the author's paper [17\ . 

4.1.4 Phase transitions 

The derivative of E{6) with respect to 6 is called the entropy, and when the entropy 
has a jump at a point 6q, it is said that a ffist-order phase transition occurs; when the 
second-order derivative has a jump, it is said that a second-order phase transition occurs, 
etc. 

But, of course, it is not the function itself that has a jump, but the leading term of its 
asymptotics as — >■ oo. In experiments, this "jump" can sometimes be actually extended 
in time, but we have agreed to consider variations in sufficiently large periods of time. 
In these periods, not only thermodynamic, but also dynamical processes can occur. We 
neglect them and consider only the time periods in which the system comes to equilibrium 
with the "energetic efficiency" taken into account [TS]. 

The author discovered the zeroth-order phase transition both in the theory of super- 
fluidity and superconductivity and in economics (a stock price break-down, a default, 
etc.), and, quite unexpectedly, it has turned out that, in view of the natural axiomatics 
(see [H]), the zeroth-order phase transition is related to quantum statistics and ther- 
modynamics. This phase transition has not been noticed by physicists, and, perhaps, it 
contradicts their ideas that the free energy can be determined up to a constant. 

4.1.5 Metastable state 

We consider a simple example of semiclassical approximation of the one-dimensional 
Schrodinger equation 

-h^y'n + U{x)yn = XnVn, Vni^) 6^2, /K 1, 

where u{x) = (x^ — lY + qx and g > is a constant. 

In classical mechanics, this equation gives the picture of "two cups with a barrier 
between them." If a particle is at the bottom of the right-hand cup with higher walls, 
then it can get into the other cup with lower walls only if the barrier disappears. 

From the viewpoint of the Gibbs postulate, as ^ 0, the particle must be in the 
deeper well. Nevertheless, it is obvious that if h <^ 1, then the particle will stay in the 
well with higher walls for a very long time. In this case, the summation in the Gibbs 
formula must be performed over the subset of eigenvalues corresponding to the well with 
higher walls whose eigenfunctions tend to zero as /i — > outside this well. Moreover, the 



temperature 6 must not be too high, so that the eigenvalues above the barrier do not play 
any role in the Gibbs formula. 

Such a state at a local minimum of the potential well u{x) is an example of a metastable 
state. 

If we consider the matrix element of the transition from the lower level corresponding 
to the well with higher walls to the lowest level Aq (corresponding to the bottom of the 
deep well), then it turns out to be exponentially small with respect to h, but any transition 
to higher levels is forbidden according to the "energetic efficiency" law (i.e., we consider 
perturbations by "friction"). 

4.1.6 Superfluidity 

N. N. Bogolyubov used Landau's ideas (see the footnote on p. 219 in [20]) to show that 
superfluidity is not a motion of the fluid, not a dynamics, but a state of the fluid such as, 
for example, ice or vapor for water. 

This state corresponds to a metastable state of the system such that any transition 
from this state to the normal state is almost forbidden: it is exponentially small as 
N ^ oo. 

N. N. Bogolyubov proved this rigorously under the assumption that the system of 
Schrodinger equations is periodic; in other words, the Schrodinger equations were consid- 
ered on a torus. The spectrum of superfluid velocities (the energy levels corresponding to 
the relevant momenta) was discrete. This readily distinguishes the state of superfluidity 
from the hydrodynamics of fluids. In the limit as the torus radius tends to infinity, the 
spectrum does not become, as usual, a continuous spectrum, but becomes an everywhere 
dense point spectrum. 

In the author's opinion, the fact that the spectrum is everywhere dense pointwise in 
the limit can be easily explained from the physical viewpoint. 

Indeed, if the system is in a state with a superfluid velocity v, then its transition 
to a state with larger velocity is forbidden by the "energy state efficiency" law, and its 
transition to a state with any lesser velocity is forbidden by the fact that any decrease 
in the velocity contradicts the notion of superfluidity. From the mathematical viewpoint, 
this means that the larger the torus radius, the less (exponentially less) is the transition 
from one state to another. 

4.1.7 Zeroth-order phase transition 

1. The author explained the spouting effect discovered in 1938 by J. Allen and H. Jones 
in which a superfluid was "flowing" through a capillary of diameter 10~^cm (in fact, this 
superfluid was at the superfluid level of a metastable state). 

The author used the two-level model to show that, at a point heated (by light) till 
the phase transition temperature 6c, the heat capacity becomes infinite, the entropy has 
a jump, and the free energy decreases to its lower value, i.e., to the point attained by 
the curve issuing from the ground state heated to the temperature 6c- This means that a 
zeroth-order phase transition occurs. In the present paper, we show that the same picture 
also appears in N. N. Bogolyubov's model of a weakly nonideal Bose gas [22]. 

This phenomenon can be easily explained if we assume that superfluidity is not a 
thermodynamic state, but the motion of a fluid. At 6 = 6c, the fluid becomes viscous and 
cannot penetrate through a thin capillary. 



But the point is that this is not a motion, but a state, and so we have a zeroth-order 
phase transition. 

The following question arises. What will happen to superfluidity not in a capillary, 
but in a rather thick pipe? Where is the zeroth-order phase transition? At ^ = 9^-, the 
fluid passes from the superfluid state to the state of the ordinary fluid with viscosity 
and begins to flow according to the usual laws of hydrodynamics. After a while, the 
motion stops, and the transition from the superfluid state to the state of a fluid at rest 
is a thermodynamic transition. All the intermediate flow is hydrodynamics, and must be 
neglected for our time scale. 

2. Let us consider the example studied in Sec. 1.5 in more detail. We slowly vary the 
constant q. We note that, in thermodynamics with the field taken into account, there are 
two more thermodynamic quantities: the field intensity and the charge. Thus, a variation 
in g is a variation in a thermodynamic variable. 

We shall show that the passage of q through the zero point results in a zeroth-order 
phase transition. Indeed, resonance occurs at g = 0: the eigenfunctions are already not 
concentrated in one of the wells and the probability (the square of each eigenfunction) 
is identically distributed over both wells. The number of eigenvalues is "doubled" and 
the Gibbs distribution has a jump. It is precisely here that one can see the role of time. 
A certain amount of time is required for half the function, decreasing outside the first 
well, to be "pumped" into the second well. But the Gibbs formula does not take this into 
account. It may happen that we will have to wait for this transition for a long time, as it 
was in the preceding example with a thick pipe. 

If q becomes negative, then all the eigenfunctions remain both in the first and in the 
second well, and the Gibbs formula is taken over all the eigenvalues, rather than over 
some of them. 

3. Finally, we consider the effect of transition into the turbulent flow for fluid helium, 
which, in fact, is very close to Landau's idea concerning energy pumping between large 
and small vortices. 

As will be shown in another paper, it is precisely the resonance that occurs between 
vortices of these two types that results in the zeroth-order phase transition, which sharply 
changes the thermodynamic parameters from the viewpoint of thermodynamics in which 
we take into account the large scales of time between transitions and, naturally, the 
averaging over these times. 

4. Since only the values of bank notes are important, while their serial numbers do 
not play any role, and the interchange of two bank notes of the same denomination is not 
an operation, bank notes obey the Bose statistics. The averaging of gains is a nonlinear 
operation, as well as addition. As was shown above, the only nonlinear "arithmetics" (the 
semiring) that satisfies the natural axioms of averaging for Bose particles (bank notes) 
has the form a © 6 = ln(e" -|- e^). 

This lead to a formula of Gibbs type. A variation in (3 = 1/T can be treated as a 
variation in the rouble purchasing power caused, for example, by printing a lot of new 
bank notes. After this, for a period of time, the balance (equilibrium) is again established. 

Ordinary financial efficiency plays the role of energetic efficiency. The zeroth-order 
phase transition is either a default or a crisis [23] . 



4.2 An exactly solvable model 

First, we consider the one- dimensional Schrodinger equation for a single particle on the 
circle 

Hijjk{x) = Eipkix), ijjk{x - L) = ipkix), (47) 

where ipk{x) is the wave function, x takes values on the circle, and if is a differential 
operator of the form 

H = e{-ind/dx), e{z) G C°°. (48) 

Here h is the Planck constant, ipk{x) = exp{ipkx), where pk = 27rhk/L, k is an 
arbitrary integer, and the corresponding eigenvalues are equal to 

Ek = e{pk). (49) 
We pass from the Hamiltonian function e{p) to a discrete function 'e{p) of the form 

e'(p) = e(nAp) for Ap(ji — < P < Ap^n + -j, (50) 



where n is an arbitrary integer and Ap is a positive constant. Then the Hamiltonian 

changes appropriately, and we denote the new Hamiltonian by H. The eigenfunctions of 
this operator coincide with ipkix), and the eigenvalues (ll9ll become 

Ek = e{Pk). (51) 
In what follows, we assume that the constant Ap takes the form 

27ihG , , 

Ap=^-, (52) 

where G is a positive integer. 

In this case, it follows from fl50l) that the set of energy levels (1511) is the set of G-fold 
degenerate energy levels 

Xn = Ecn = e{pGn)- (53) 

The Schrodinger equation for noninteracting particles has the form 

Hn'^{,xi, . . . , xtv) = ^^(xi, . . . , xtv), (54) 

where '^{xi, . . . ,xn) is a symmetric function of the variables xi, . . . , xat (bosons). The 
Hamiltonian is given by the formula 

H^ = Y,Hj, (55) 
i=i 

where Hj is the Hamiltonian of the particle with number j, which has the form 

The complete orthonormal system of symmetric eigenfunctions of the Hamiltonian (!55|) 
has the form 

<I/{n}{xi,...,xn) = Symm^^^ ^^^ilj^N}{xi, . . . ,xn), (57) 

Ifi 



where Symm^^ is the symmetrization operator over the variables xi, . . . ,xn, {N} is 
the set of nonnegative integers Nk, A; G Z, satisfying the condition 

oo 

J2 Nk = N, (58) 

A;=— oo 

and the function ip^Ny^xi, . . . ,xn) is equal to 



N 



'iIj{n}{xi,...,xn) = Y\_^kAxs)- (59) 

s=l 

Here the indices ki, . . . ,kN are expressed in terms of the set {A^} using the conditions 

N 

ks < ks+i for all 1 < s < iV - 1, ^ 6kk, = Nk for all k e Z, (60) 

s=l 

and 6kk' is the Kronecker symbol. The eigenvalues of the Hamiltonian (|55l) are 

oo 

E{{N}) = J2 ^k^k- (61) 

A;=— oo 

We consider interparticle interactions of the following type. We assume that the particles 
interact in pairs and the interaction operator for particles with numbers j and k has the 
form 

% = -^W{H,-Hk), (62) 
where \^ > is the interaction parameter and the function W{E,) is given by the formula 

. X fl for If I < D, 

^ lo to. !«! > D. <«^' 

Here > is the parameter of the width of interaction with respect to energy. The 
operator (jHS]) corresponds to the interaction under which particles in a pair attract each 
other and radiate the energy quantum —V/N if the difference between their energies is 
less than D and do not interact at all if the difference between their energies is larger 
than D. The Hamiltonian of the system of particles with interaction (163|1 has the form 

N ^ N N 

i=l j=l k=j+l 

In view of fl62|) . the sums in flM|) commute, and hence the set of eigenf unctions of the 
Hamiltonian ( 164|) coincides with ( 1571) . It also follows from (162|) that the corresponding 
eigenvalues are equal to 

oo oo oo 

E({N}) = EkNk -^Y. H "^(^k - Ei)(NkNi - 4iNk). (65) 

k=— oo k=— oo 1=— oo 

In what follows, we assume that the interaction width is sufficiently small, satisfying the 
condition 

D < min |A„ — Am|. (66) 
17 



The set of energy levels Ek, k E Z, coincides with the G-fold degenerate set of levels fl53l) . 
Hence, by ( l66i) . the set of energy levels ( l65l) of the system of particles under study can 
be written as 

OO -j^ oo 

E({iV}) = E - ^ E - 1)' (67) 

n=—oo n=—OQ 

where the level E{N}) has the multiplicity 

r({;v}) ^ n (68) 

nil (G-l)!iVJ 

Here {A^} denotes the set of nonnegative integers A'^, n G Z, satisfying the condition 

OO 

J2 ^n = N. (69) 

n=— OO 

Let us consider the partition function for the system of bosons with Hamilto- 
nian (16^ . Since the energy levels and their multiplicities are given by formulas (!67|) 
and fl68p . respectively, the partition function at temperature 6 takes the form 

Z{e,N) = Y,mN})exp{{-E{N})e). (70) 

{TV} 

Here the summation is performed over all sets {A^} with condition fl69|) taken into account. 

In what follows, we assume that G depends on A^ and the following condition is 
satisfied: 

G 

lim — = ^ > 0. (71) 
By F{{N}, 9) we denote a function of the form 

F({iV}, 6) = E({N}) - e ln(r({N})) , (72) 

and by {A^°} the set of nonnegative numbers A^°, n G Z, for which the function (172|) is 
minimal under condition (l69|l . 

Now, we consider the problem of finding the minimal value of the function fl72p under 
condition (1691) . In the limit as A^ — >^ cxd and under condition 0711) . the point of minimum 
has the form 

Ar„(0) = Ar(m„(0) + o(l)), (73) 
where m„(6'), n G Z, is determined by the system of equations 

(TTl \ 
— = fi{9), neZ, (74) 
g + runj 

and fi{9) can be found from the equation 

OO 



J2 ^n = l- (75) 



n=— OO 
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The substitution of fl73l) into fl72l) and then the use of the asymptotic Stirhng formula 
give the following relation for the specific free energy: 



fie) ^ hm fie, N) = hm :WV^ ^^^^ 

= E (^Xnmr,-^m^}j + (^em^\n(^^^ (77) 

where, for brevity, we omit the argument e of rrinie). We introduce the notation 

= A„ - Fm„. (78) 
Remark 3. Remark In the notation fl78l) . the system of Eqs. fl7^ and fl75l) takes the form 

u^nie) =Xn- V ^ . ^ , n G Z, (79) 

exp((u;„ - fi)/e) - 1 

oo 

exp((.;„-/.)/^^)-l = 

The system of Eqs. (17^ and (ISUl) exactly coincides with the temperature Hartree equa- 
tions [23] for the system of bosons with Hamiltonian (1641) . 

We shall study the solutions of Eqs. (!74|) and (!75|) . For = 0, the system has many 
solutions, which we number by the integer /: 

mll^ = Sin, nJeZ. (81) 

Among all the numbers /, we choose those that satisfy the condition 

Uni = Xn-Xi + V>0 for all n I. (82) 

For these numbers, there exist solutions of Eqs. fl74p and (1751) converging to (1811) as ^ ^ 0. 
The asymptotics of these solutions as 6* ^ has the form 



m^^^ gexpi-^j \fn ^ I, 1 - mf^ ~ ^ ^ exp ( j . (83) 

Thus, depending on the spectrum A„, n G Z, for sufficiently small values of the tempera- 
ture e, the system of Eqs. fl7i|) and fl75|) has many solutions. These solutions, along with 
the point of global minimum, also contain points of local minimum of the function (l72l) . 
The values of the function ( 172|) at the points of local minimum are equal to the free energy 
of metastable states. We consider the function ( 1721) at 6^ = 0. In this case, it coincides 
with the energy spectrum of system (1671) . We consider the energy of the system for the 
case in which almost all particles are at the energy level A;, which means that following 
conditions holds: 

iV„ < Vn ^ /. (84) 

Deriving Ni from relation ( |69l) and substituting the result into (|67l) . we see that, in view 
of (|8^ . the energy spectrum of the system in this domain has the form 

E(N) ^ AiN - ^ + V(A„ - Ai + V)N„. (85) 

IP 



To the Hamiltonian flMl) . there correspond the Hartree equation and the appropriate 
system of variational equations. To each / G Z, there corresponds a solution of the Hartree 
equation, and this solution describes the state 

iV« = N5ni. (86) 

Moreover, the eigenvalues of the system of variational equations corresponding to this 
solution of the Hartree equation coincide with 

T^ni = \n-k + V, n^l. (87) 

In [25l [26] . it was shown that if the eigenvalues of the system of variational equations 
for the solution of the Hartree equation are real and nonnegative, then such a solution 
corresponds to the ground state or to a metastable state of the system. This means that 
to all / for which condition fl82|) holds at ^ = 0, there correspond metastable states of the 
system of bosons. As in the case of zero temperature, to solutions of Eqs. fITil) and fl75|) 
for 6* 7^ 0, there corresponds a temperature metastable state if the point ( |73|) is a point of 
local minimum. Now, we note that, for very large temperatures, the system of Eqs. ( TMl) 
and ( 1751) has only one solution corresponding to the global minimum of the function ( 1721) . 
The asymptotics of this solution as ^ — > cx) has the form 

nm{0) ~ g 



The uniqueness of the solution at large temperatures means that all metastable states 
disappear with increasing temperature. The temperature at which a metastable state 
disappears is critical for this state. 

Let us consider the behavior of the entropy and the heat capacity of metastable states 
when we approach the critical temperature. We consider the metastable state to which 
there corresponds a solution of Eqs. (17^ and (175|) converging to (IHTj) as 6* — > for some / for 
which ( l82l) holds. In what follows, we assume that the function e{p) satisfies the condition 
e(0) < e{p) for p ^ 0. Then, according to ( l53l) . we have Aq < A; for / 7^ 0. Hence the 
solution of Eqs. ( 174|) and ( 1751) converging to ( IHTl) as 6* — >■ for / = corresponds to the 
temperature ground state of the system of N bosons. Moreover, condition (l82l) becomes 
equivalent to the condition A^ — Aq < V. 

We assume that this inequality holds for 1^0. The condition that the corresponding 
solution rrin^ (6) of Eqs. (I74p and (I75p determines a point of local minimum of the func- 
tion (1721) under condition (1691) by formula (I73p can be written as the following system of 
inequalities: 



wi^\e){g + m'i\e)) 
- -v + ^, a. < 0, -E^r^ < 1 



We note that these inequalities hold for ([HSD as — 0. Inequalities (15^ follow from 
the condition that the second variation of the function (1721) is positive; the variation is 
calculated under condition (l69l) . 

The metastable state disappears at the temperature at which the last inequality in (189]) 
becomes an equality. We denote this critical temperature by 6c \ It follows from (189|) and 

9n 



Eqs. (17^ and (1751) for 6 < 6c^ that mn\0) is an increasing function of the variable 6 for 
n ^ I and mf\e) is a decreasing function of 6'. Moreover, we see that 

mf\e) > mi^{e) > mJV) (90) 

if A.„ < A„/ and n,n' I. 

From (1771) we obtain the following expression for the specific entropy of a metastable 
state in the limit as ^ oo: 

^^'\(^)= E U9 + mll^)\M + '^]-rn^hnl^]\, (91) 

n=— oo ^ ^ 9 / \ 9 / / 

where, for brevity, we omit the argument 6 of mn\0)- Differentiating (PTl) . we obtain 
The last inequality follows from the properties of mn\6). Since the last inequality in (l89l) 



becomes an equality at the critical temperature, we can show that, as 6 Oc'^ — 0, the 
solutions of Eqs. (17il) and ( 1751) behave as follows: 

r'W / 

mll\e) - rr^ie^) ^ -^jv-nr^\J e^P -e g z, (93) 

On (6'c ) 

where C*^'^ is a negative number. The substitution of (193|) into (19T|) shows that the 
derivative of the specific entropy with respect to the temperature (this quantity is equal to 
the heat capacity divided by the temperature) tends to infinity as 6* — oP — according 



to the law 1/\J6P — 6. This means that the projection of the Lagrangian manifold 
corresponding to the metastable state on the ^-axis is not uniquely determined in a 
neighborhood of the critical temperature [HI 127]. The derivative of the temperature 
with respect to the entropy vanishes as the critical temperature is approached. Therefore, 
as well as in view of (l92l) . as was already pointed out, the Lagrangian manifold is uniquely 
projected on the s-axis. We note that it follows from the properties of rrSn(6) that the 
following conditions hold for 6 < oP: 

m«(^)<m«(^) Vn^O,/, 



91 



Inequalities flMj) mean that the potential barrier between the energy levels A/ and Aq 
is less than the energy barrier between A/ and A„ for n ^ 0, /. This means that, as 
the critical temperature is attained, the system of bosons under study which is in the 
metastable state with number / changes its state in a jump and passes to the temperature 
ground state. This is a zeroth-order phase transition, since not only the entropy and heat 
capacity, but also the free energy have jumps. 

If the quantity 

min |A„ — Ao| =5 (95) 

is small, then the difference between the free energies in the zeroth-order phase transition 
from the lowest metastable state to the ground state is also small, and the heat capacity 
in this transition has a singularity. The asymptotics of the partition function near the 
critical point is given by the asymptotics of the canonical operator in a neighborhood of 
a focal point. This asymptotics has the form of an Airy type function of an imaginary 
argument and has a singularity as ^ oo. The parameters 5 <^ 1, ^ 1, and L ^ 1 
can be chosen so that the form of the phase transition point coincides exactly with the 
A-point [2H]-[33]. 

If adequate measures increasing the Kolmogorov complexity of the Hartley entropy 
are not taken, the financial system still does not necessarily result in a zero-order phase 
transition (social explosion). 

Now, let us return to the problem about deposits in a pyramid and banks, i.e., to 
the assumptions that the bills of the same value with different numbers are identical. 
The phase transition related to the disappearance of the condensate has the following 
meaning in this case. As /? (price) decreases, the actions cease to be bought and sold 
at some moment /3o, i.e., nobody trades them aX. (3 < (3q, although, seemingly, it is more 
profitable to sell them at any price, and therefore, somebody can speculate in resale, 
thereby reducing the price of the actions to zero. Nevertheless, this paradoxical fact is 
observed in practice and mentioned in the literature. 

Then the corresponding shares will be involved in a complicated barter exchange, 
which makes for an increase in the entropy. The same happens to other equities, in par- 
ticular, to currencies provided that the inflation rate is sufficiently high. All in all, this 
ensures that the entropy and complexity increase and can prevent zero-order phase transi- 
tion, which occurs only at some threshold temperature and for some nonlinear interaction 
swinging the equilibrium system. 

Now we consider the following model. Assume that the number of people is = ni+n2. 
The culture (velocity, energy) level Ai typical of rii is twice less than the level A2 typical 
of n2. The quantities Ai = 1, A2 = 2, and 1 < 7 < 2 determine the quadratic interaction 

where ni of people have A^i units of money and ^2 of people have A"2 units of money, 
and Ni N2. For a sufficiently large M/Nq (an analog of temperature), there is a phase 
transition to the state with Ni = N and A"2 = with an enormous outburst of the kinetic 
energy (in physics, this is the Allen- Jones fountain effect [H fT7]). 
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